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Abstract: We study the light scattering of homogenous radially-anisotropic 
spherical particles. It is shown that radial anisotropy can be employed to 
tune effectively the electric resonances, and thus enable flexible overlapping 
of electric and magnetic dipoles of various numbers, which leads to unidi¬ 
rectional forward super-scattering at different spectral positions. We further 
reveal that through adjusting the radial anisotropy parameters, electric 
and magnetic resonances of higher orders can be also made overlapped, 
thus further collimating the forward scattering lobes. The ultra-directional 
super-scattering we have obtained with individual homogenous radially 
anisotropic spherical particles may shed new light to the design of compact 
and efficient nanoantennas, which may And various applications in solar 
cells, bio-sensing and many other antenna based researches. 
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1. Introduction 

Recently stimulated by the demonstrations of optically-induced magnetic responses of high 
primitivity dielectric particles GM, the original proposal by Kerker la and the concept of 
Huygens source in the antenna theory both of which are related to the manipulation of 

particle scattering patterns, have attracted renewed surging interest 181 42^ . The central precon¬ 
dition for the realization of Kerker’s proposal and Huygens source is the efficient excitation 
of magnetic resonances and their interferences with the electric counterparts EM). Such in¬ 
terferences between electric and magnetic resonances can be employed for efficient shaping of 
scattering patterns for nanoparticles in the far-held Il8l- l26l . which thus might paly a critical role 
in many applications for biological sensing, compact nanoantenna design, photovoltaic devices 
and so on. 

Among all kinds of manipulations of the scattering patterns based on optically-induced mag¬ 
netic responses, simultaneous forward scattering enhancement (forward super-scattering) and 
backward scattering suppression based on overlapping electric dipoles (EDs) and magnetic 
dipoles (MDs), which satishes the hrst Kerker’s condition BET), is one of the most attractive 
topics UMBEHSl- Though for homogenous dielectric spheres, both EDs and MDs can be 
excited and actually unidirectional forward scattering with suppressed backward scattering has 
been experimentally demonstrated iminiiiii. The problem is that the corresponding func¬ 
tioning wavelength is far from the resonant position, due to the fact that the EDs and MDs 
are spectrally separated, resulting in reduced overall scattering. To obtain simultaneous unidi¬ 
rectional scattering in the super-scattering regime ll28H30l . with enhanced forward scattering 
and sufficiently suppressed backward scattering, basically EDs and MDs have to be spectrally 
overlapped. This has recently been realised in core-shell plasmonic nanoparticles ll^ 1 01119113IL 
dielectric nanodisks ll20l or in spheroidal dielectric nanoparticles ll24ll . 

Here in this paper, we demonstrate an alternative way based on radial anisotropy to spectrally 
overlap electric and magnetic resonances, which subsequently leads to unprecedented scatte¬ 
ring patterns for individual homogenous particles. It is shown that in spherical dielectric parti¬ 
cles, the resonant positions of electric multipoles can be effectively tuned by radial anisotropy, 
thus enabling flexible overlapping of electric and magnetic resonances of various orders and 
numbers. Beyond the conventional scattering manipulation to obtain unidirectional forward 
super-scattering with negligible backward scattering based on overlapped EDs and MDs, we 
demonstrate that the directionality of the forward scattering can be further improved. This has 
been realized by overlapping electric and magnetic resonances of higher orders through opti¬ 
mizing radial anisotropy parameters, which result in ultra-directional forward super-scattering 
of individual homogenous particles. The mechanism we have revealed of radial-anisotropy in¬ 
duced ultra-directional forward super-scattering can be certainly extended to particles of other 
kinds of anisotropy Il32l434l and shapes, and it is expected that such mechanism can shed new 
light to many particle scattering based applications, especially those which require flexible scat¬ 
tering pattern shaping and control. 






Fig. 1. (a) Scattering of an incident plane wave by a radially anisotropic spherical parti¬ 
cle of radius R, radial permittivity Er and and transverse permittivity £t = . The radial 

anisotropy parameter is defined as T] = e, /Er. The plane wave is polarized (in terms of elec¬ 
tric field) along x direction and is propagating along z direction, (b) Normalised scattering 
cross sections for ED [Nsca(‘j>i), red curves] and MD [Nsca(^i)i blue curve] of different 
anisotropy parameters. The resonance positions of EDs and MDs are indicated hy Eii,i 2 
and Mil 12 respectively, with the corresponding near-held distributions at those points [on 
the x — z plane of y = 0 in terms of both out of plane H^, (color-plots) and in plane E (vector- 
plots)] shown in the bottom row. The dashed black curves indicate the boundaries of the 
particles, (c) The dependence of resonance positions [solutions of Eq. ([6ll] of EDs (eiLi2) 
and MDs (mu 12 ) on the anisotropy parameter r\. 







































































2. Theory of plane wave scattering by radially anisotropic spherical particles 

As is illustrated in Fig. [TJa), we study the scattering of plane waves by radially anisotropic 
spherical particles (with radius R). The effectively permittivities along radial and transverse di¬ 
rections are £,■ and e, = respectively, with the anisotropy parameter defined as Tj = e, /Cr. The 
plane wave is polarized along x direction and propagating along z direction. The scattering can 
be solved analytically using modified Mie theory Il33l43^ and the total normalised scattering 
cross section Nsca [normalized by where k is the angular wave number in the background 
(it is vacuum in this study)] is: 


A^sca= + (1) 

n=\ 

where and are Mie scattering coefficients and for radially anisotropic particles can be 
expressed as 1(3^143^ : 


_ m\i/fi{ma)\i/'„{a) - v/„(a)v4(wa) 
^ _ v/„(wa)v/'(a) -mv/n(a)v/'(OTa) 

where 


n = ^n{n + l)ri + ^-^, 


( 2 ) 

(3) 


(4) 


and a = kR\ \j/ and ^ are Riccati-Bessel functions ll35l and the accompanying primes indicate 
their differentiation with respect to the entire argument. Considering that a„ and correspond 
to n-th order electric and magnetic multipoles respectively {e.g., ai and bi correspond to ED 
and MD respectively), then the Asca contributed from them can be expressed respectively as: 


Asca(a«) = (2u +l)|a„p, Asca(^«) = (2n-f l)|h„|^. (5) 

Since |a,i|, |h„| < 1 03311351 . the maximum of a single n-th order multipole (either electric 
or magnetic) would be 2n + 1, where super-scattering is defined as the case when the total Nsca 
is large than this number due to the coexistence to multiple multipoles li29l . The resonance 
positions of the electric and magnetic multipoles of the n-th order are defined as the solutions 
(in terms of a) of the equations below respectively: 

\an (Cnj) I = 1, \b„ (mnj) I = 1, (6) 

where j is the number of the solution [counted from smaller to larger a (from smaller to larger 
particle radius, or from larger to smaller wavelength)] and thus corresponds to the number of 
a specific multipole resonance llJ7l439l . For example, a = ei 2 is the resonance position of the 
second ED, and we would like to note that most of the previous work only focuses on the 
case of 7 = 1 (resonances supported at largest wavelengths) and neglect electric and magnetic 
multipoles of higher numbers (supported at smaller wavelengths) l[TI44]l9] 424ll . 

According to Eq. (|2|i-Eq. Q, magnetic multipoles (characterized by b„) will be dependent 
on transverse permittivity £( = only and will not be affected by the radial anisotropy if 
m is fixed. In contrast, electric multipoles (characterized by a„) is dependent on the radial 
anisotropy [through «; see Eq. ®] and thus can be tuned. Those different features of elec¬ 
tric and magnetic multipoles can be explained as follows: magnetic resonances correspond to 





transverse electric-field (TE, electric field is perpendicular to the radial direction) distributions 
as shown in Figs. [T]-Mii 12 , which means that the change of the radial permittivity would not 
at all affect the magnetic resonances; Electric resonances however, correspond to transverse 
magnetic-field (TM, magnetic field is perpendicular to the radial direction and there are both 
radial and transverse electric field components) distributions, which indicates that the change 
of the permittivity in the radial direction would influence the electric resonances. Except the 
explicit expressions of a„ , other expressions for radially anisotropic spherical particles (h„ and 
those expressed by and h„), including extinction and absorption scattering cross sections, 
angular scattering intensities, and the scattering matrix will be the same as those for isotropic 
spherical particles 03313511 . As a result, all the scattering pattern manipulation approaches, espe¬ 
cially those based on interferences of electric and magnetic multipoles Il8l 424ll . can be directly 
applied for radially anisotropic particles. 

3. Unidirectional super-scattering of radially anisotropic particles with overlapped EDs 
and MDs 

In Fig.lHb) we show the Asca of both EDs and MDs. Throughout the paper without losing the 
generality, we fix m = 3 .5 and change the radial permittivity to tune anisotropy parameter rj , and 
consequently the resonance positions of MDs will be fixed. As described in the section above, 
there are multiple EDs indicated by Eq and MDs indicated by Mq (here we have only shown 
mode number J up to 2 only, and EDs and MDs of higher numbers can be supported at larger a). 
The corresponding near-held distributions are shown (on the x-z plane with y = 0) in the bottom 
row of Fig. [T] where the color-plots correspond to the perpendicular magnetic helds along y (Hy) 
and the vector-plots correspond to the electric helds on the x-z plane. It is clear that for MDs 
there are effectively current (denoted by electric held) circulations with strong magnetic held 
distribution at the center of the circulation (Figs. [BMii, 12 ); and for EDs there are effectively 
linearly polarized currents with negligible magnetic held distribution at the center of the particle 
(Figs.IBEii, 12 ). As is also indicated by the near-held distributions, the mode number J basically 
represents the number of anti-nodes of the helds from the centre to the boundary of the particle, 
which is highly related to the quantum-mechanical shape resonances ISHISI. We note here 
that for resonances of the same order n, different mode order J can be only recognized from 
the near-held distributions. In the far-held however, mode of the same order but different mode 
numbers are indistinguishable. 

Nextly, we study the effect of radial anisotropy on the resonant positions and show the results 
in Fig.[Tlb). As has already been described in the section above, with hxed m MD will not be 
affected by the radial anisotropy, while in contrast ED will red and blue shifted with decreasing 
and increasing 77 respectively [Fig. [Ub)]. This means that basically the relative resonant posi¬ 
tions of EDs and MDs can be tuned through change 77 , and moreover we can overlap hexibly 
EDs and MDs even for different mode numbers. 

In Fig. me) we show the 77 dependence of central resonant positions of EDs (eq) and MDs 
(mq) for mode number up to 2 only. Obviously there are three crossing points, where ED 
and MD are overlapped. This can be confirmed by their corresponding normalized scattering 
cross section spectra shown in Figs. [2ja)-l2tc) respectively, where the total cross section has 
also been shown by black curves. For all three cases, ED and MD are perfectly resonantly 
overlapped, and the difference for the three cases is the number j of the modes overlapped. 
Previous studies |I9HII][I1[I51|IM01|23]|241 focuses only on the case of j = I (resonance at 
the largest wavelength) while here we have demonstrated that EDs and MDs of higher mode 
numbers (supported at smaller wavelengths) can be also mode overlapped, providing more 
freedom for scattering pattern engineering based on interferences between electric and magnetic 
multipoles. 
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Fig. 2. Normalised scattering cross section spectra for overlapped EDs (red curves) and 
MDs (blue curves): (a) y = 1 for both ED and MD with rj = 0.103; (b) J = 2 for both 
ED and MD with rj = 0.179; (c) j = 1 for ED and j = 2 for MD with T] = 5.252. Eor 
all three cases the total normalised cross section spectra have been provided (black curves) 
with the overlapping resonant positions of a = 0.861, 1.821, 1.821 indicated by black dots 
respectively. The corresponding two-dimensional (red curves: scattering on the x — z planes; 
blue crosses: scattering on the y — z planes) and three-dimensional (a part is cut off for better 
visibility) scattering patterns are shown in the middle and right columns respectively. 


The corresponding scattering patterns for the three cases at the overlapping resonant positions 
[indicated by black dots in Figs. |2ta)-l2jc)] are shown in the middle column [two-dimensional 
scattering patterns on the x — z (red curves) and y — z (blue crosses) planes] and right column 
(three-dimensional scattering patters with a part cut off for better visibility) of Fig.|2] For the 
case of Tj = 0.103 [Fig.|2l[a)] at the resonant position (a = 0.861), since the total Nsca (which 
equals approximately to 6 ) is about twice of that of the (single-dipole limit, which is 3), mean¬ 
ing that scattering contributions from other multipoles are negligible. This lead to the more or 
less ideal case of unidirectional forward super-scattering (with null backward scattering and the 
scattering pattern is azimuthally symmetric) from overlapped ED and MD |l5]|3[T0l, as demon¬ 
strated inFigs.|2a-l) and|2ta-2). For the other two cases of 7] =0.179 [Fig.|2l[b)] and 77 =5.252 
[Fig.EJc)], since the contributions from other multipoles are not negligible (total Nscn notably 
larger than 6 at the overlapping resonant position of a = 1.821), the scattering is not ideally 
azimuthally symmetric and moreover there are some minor backward scattering [see Figs.^b- 
l);|2tb-2),|2jc-l) and|2jc-2)]. Overall however, for all three cases the backward scattering has 
been significantly suppressed and forward scattering has been enhanced, justifying our claim 
of unidirectional forward super-scattering of radially anisotropic particles. 












4. Ultra-directional super-scattering of radially anisotropic particles with overlapped 
electric and magnetic multipoles of higher orders 

It is recently shown that through overlapping electric and magnetic multipoles of higher orders, 
ultra-directional forward scattering pattern could be obtained II23I40I41II . Those demonstrations 
either rely on specific incident dipole sources IMP or on composite core-shell structures ll 2 ^ . 
where the condition of super-scattering cannot always be guaranteed. Here we try an alternative 
way to employ radial anisotropy to overlap electric and magnetic multipoles of higher orders, 
thus obtaining ultra-directional forward super-scattering of homogenous particles with plane 
wave incidence. 

Without losing generality here we confine the discussion to electric and magnetic multipoles 
of mode number 7 = 1. First of all to directly evaluate the radial anisotropy parameter rj required 
to overlap electric and magnetic multipoles of different orders, we have extended the simplified 
expressions for the central resonance positions of EDs and MDs in isotropic particles llJTl to 
radially anisotropic particles as follows; 

_ Z 1 3 1 

m 7}l^m (m^ — 1)'/2 10 3 — 1)^/^ ’ 

, s X AiX^/3 _ 1 3 1 AiX-^/3 

m l}l^m —1)^^ 10 2^l^m 3 2F3(^;7j2 _ jp/2 ’ 

where Ai w 2.338 is the first root of the Airy function of Ai(-x), = « -|- 1/2 and x=n+l/2. 

The overlapping of electric and magnetic multipoles of the n — th order requires that: 


(7) 

( 8 ) 


^nl(t7nl) ^nl —0, 


(9) 


where tlni is the solution and thus exactly the anisotropy parameter required. 

In Fig. 13 a) we show up to mode order of n = 10, where both approximated [dashed 
blue curve, Eq. 0 -Eq. 0] and exact results [red crosses, calculated based on Mie theory; see 
Eq. 0-see Eq. (|6]l] are provided. In the inset of Fig. [3a) we show the corresponding Nsca 
spectra for electric and magnetic multipoles of order n = 6, confirming the ideal multipole 
overlapping. The corresponding two-dimensional scattering patters for overlapped quadrupoles 
(77 = 0.3747), hexapoles (77 = 0.5254), and multipoles of 77 = 10 (77 = 0.8269) are shown in 
Figs.|3b)-|3d) respectively. As has been shown, the backward scattering has been sufficiently 
suppressed [the minor backward scattering and azimuthal-asymmetry originate from the extra 
contributions from other multipoles, which is similar to that shown in Fig.lH and the forward 
scattering has been more and more collimated with increasing overlapping multipole orders. To 
characterize quantitatively the forward scattering directionality, we define the angular beam- 
width j3 of the main forward scattering lobe, as shown in Fig.|3b), which correspond to the full 
width at half maximum of the scattering intensity. We show the dependence of j3 on mode order 
in Fig. [3e) and it is clear that as mode order increases, the forward scattering directionality 
can be further improved. This confirms the achievement of ultra-directionality of the scattering 
from individual homogenous anisotropic particles. This has been further visualized by the insets 
shown in Fig. |3e)> which show three dimensional scattering patterns of overlapped multipoles 
of « = 2,3,10. 

5. Conclusions and discussions 

To summarise, we study the scattering patterns of radially anisotropic homogenous spherical 
particles with incident plane waves. We demonstrate that radial anisotropy can be employed 
to tune effectively the central resonance positions of EDs, and thus alow flexible overlapping 
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Fig. 3. (a) The anisotropy parameter required (T]ni) to overlap electric and magnetic mul¬ 
tipoles of mode order up to n = 10 with fixed mode number 7 = 1. Both approximated 
[dashed blue curve, Eq. 0-Eq. 0] and exact results [red crosses, calculated based on Mie 
theory; see Eq. 0-Eq. 0] are provided. Inset: normalised scattering cross section spectra 
for overlapped electric (red curve) and magnetic (blue curve) multipoles of n = 6. Two- 
dimensional scattering patterns are shown (red curves: scattering on the x — z planes; blue 
crosses: scattering on the y — z planes) for overlapped electric and magnetic quadrupoles 
(b), hexapoles (c) and multipoles of u = 10 (d), where the corresponding anisotropy param¬ 
eter Tj = 0.3747, 0.5254,0.8269 respectively. As is shown in (b), the angular beam-width 
P corresponds to the full width at half maximum of the scattering intensity of the main 
scattering lobe, (e) The dependence of (5 on the order of the overlapped multipoles and the 
insets show the three-dimensional scattering patterns of overlapped electric and magnetic 
quadrupoles, hexapoles and multipoles of n = 10. 


















of EDs and MDs of various mode numbers, which results in unidirectional forward super¬ 
scattering patterns with significantly suppressed backward scattering. We further reveal that the 
forward scattering lobe can be better collimated by optimizing the radial anisotropy parameters 
to overlap electric and magnetic multipole of higher orders, without compromising the feature 
of significant backward scattering suppression. 

We note that for the overlapping of multipoles of higher orders, we confine our discussions to 
mode number j = \ and such investigations can certainty be extended to multipoles of higher 
mode numbers, as we have done for dipoles of different mode numbers shown in Fig. [Hand 
Fig . m Also here in this paper we have discussed only radial anisotropy (for natural materials, 
some anisotropy parameters we have adopted in this paper are not quite practical, which how¬ 
ever could probably be realised in artificial metamaterials where a wide range of exotic effective 
permittivities and permeabilities can be obtained; see e.g.. Ref. iHSl l and thus it is natural to 
expect the employment of anisotropy of other types II32II33II (including magnetic anisotropy) 
would render more flexibilities for resonance tuning and thus scattering pattern manipulations. 
Moreover, when lossy anisotropic materials in included, the effects of loss would be twofold; 
(1) the resonance positions of both electric and magnetic multipoles would shift 1^ , and con¬ 
sequently the anisotropy parameters required to overlap them would be different from those of 
lossless structures; (2) material loss would change the magnitudes of both electric and magnetic 
multipoles, which consequently would not necessarily be equal to each other at the overlapping 
point, and thus the backward scattering cannot be ideally suppressed ll9l [T0ll2^ . The scattering 
shaping based on anisotropic nanoparticles has been discussed in Ref. Il34ll . which however is 
contrastingly different what we have presented in this work: (1) They confine their discussions 
to the quasi-static limit, while our discussions are more general, and can be applied to parti¬ 
cles of sizes and functioning spectral regimes that are beyond the quasi-static regime; (2) They 
confine their discussions to overlapping electric and magnetic dipoles only, while we have dis¬ 
cussed overlapping resonances of different orders and numbers, which can further collimate the 
forward scattering pattern without compromising the feature of backward scattering suppres¬ 
sion; (3) Their demonstrations are highly dependent on the magnetic responses and/or magnetic 
anisotropy, while our demonstrations are based on electric responses and/or electric anisotropy 
only and thus relatively more practical. 

Moreover, we should keep in mind that the ultra-directional forward scattering we have ob¬ 
tained is based on individual particles, and the directionality can be further improved by ar¬ 
ranging the anisotropic particles in arrays imnoi. Here we focus on the forward scattering en¬ 
hancement with backward scattering suppression (namely the first Kerker’s condition ll5l lT7ll l. 
and radially anisotropy can certainly applied to suppress the forward scattering with enhanced 
backward scattering (namely the second Kerker’s condition ||5]|23), or to control the scattering 
intensities at other scattering angles |[T^ . The mechanism we have revealed in this paper is 
not constrained to spherical homogenous particles and can certainly be extended to particles 
of other shapes and/or of many layers, which is quite promising for various particle scattering 
related applications in the fields of nanoantennas, solar cells, bio-sensing and so on. 
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